Large D — 2 Theory of Superconducting Fluctuations in a Magnetic Field and Its 

Application to Iron-Pnictides 
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A Ginzburg-Landau approach to fluctuations of a layered superconductor in a magnetic field is 
used to show that the interlayer coupling can be incorporated within an interacting self-consistent 
theory of a single layer, in the limit of a large number of neighboring layers. The theory exhibits two 
phase transitions - a vortex liquid-to-solid transition is followed by a Bose-Einstein condensation 
into the Abrikosov lattice - illustrating the essential role of interlayer coupling. Using this theory, 
explicit expressions for magnetization, specific heat, and fluctuation conductivity are derived. We 
compare our results with recent experimental data on the iron-pnictide superconductors. 



The discovery of high-temperature superconductivity 
in iron-pnictides [3, Q has led to a renewed interest in the 
physics of layered compounds and the role of supercon- 
ducting fluctuations. In older high-T c superconducting 
cuprates, due in large part to their extreme anisotropy, 
the fluctuations have taken center stage, particularly in a 
magnetic field 0]. At present, a rather good understand- 
ing of such fiuctations is available in two-dimensional 
(2D) and three-dimensional (3D) systems. However, the 
intermediate regime, where the interlayer coupling is too 
weak to be ignored and yet not strong enough to render 
the system fully 3D, remains an important challenge. Al- 
though most theoretical models of pnictides so far have 
focused on the 2D nature of these materials , exper- 
imental evidence frequently suggests a pronounced quasi 
3D behavior 0, U3], especially within the so-called 122 



family 11]. Thus, the iron-pnictides apparently belong 



to this in-between regime. 

In this Letter, we introduce a theoretical approach that 
allows for an explicit approximate solution to the prob- 
lem of superconducting fluctuations in this challenging 
intermediate situation. First, we show that the Joseph- 
son coupling between superconducting layers in a mag- 
netic field can be recast as a contribution to the effective 
"on-site" Ginzburg-Landau (GL) free energy of a single 
layer, in the limit of a large number of neighboring lay- 
ers. The system is thus described by an effective 2D GL 
theory, which - for practical purposes - can be treated ex- 
actly, by solving a set of non-linear, self-consistent equa- 
tions, in combination with a solution for the purely 2D 
case [I3 - H5j . Second, we show that this theory - unlike 
the 2D one - possesses two phase transitions, reflecting 
the crucial role of Josephson coupling. Finally, we apply 
our theory to study fluctuation effects around the upper 
critical field H C 2(T) and compare the results to recent 
experimental data on the iron-pnictide superconductors. 

We consider a general Josephson-coupled layered sys- 
tem, with an individual layer described by the GL model. 
The partition function is 



-So(i)-E, m s int (i,j) 



where ipi G LLL is the fluctuating GL order parameter 
in the ith layer; LLL denotes the lowest Landau level 
for charge 2e; and j is summed over nearest neighbors of 
layer i. The corresponding action is 



So(i) 
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where s is the distance between layers, and a = ao(t — 
t c2 (h)), where t = T/T c (0) and h = H/H c2 (0) are the di- 
mensionless temperature and magnetic field, respectively. 
The interlayer portion of the GL-LLL action ([I} is [Tjj 
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E 
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d 2 r| [^(r)^-(r)+^(r)^(r)l3) 



The goal now is to integrate out the Josephson-coupled 
portion and obtain a partition function for the 0th layer 
that is entirely "local," i.e. defined on a single layer. As 
a first step, we assume that this can be done for the lay- 
ers (denoted by j) that are adjacent to the 0th layer, i.e. 
that all couplings Si n t(j, j + c), where a denotes all layers 
neighboring layer j except for the 0th layer, can be inte- 
grated over, giving a correction to the "on-site" action, 
so that So(j) — > S' (j). (When the number of layers j is 
very large they decouple from each other, and we are left 
with a Bethe lattice, where each lattice "site" is actually 
a 2D superconducting layer and the coordination number 
of the lattice is d. This is different from Ref. [17], where 
each site is a 0D quantum cluster.) We obtain 



Z(0) = I D(V>o,V>o)e 

d 



-S o (0) 
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where Zo(j) = Z(j)\s int =o- Expanding the interlayer 
term in (j4|), and noting that only even terms in the ex- 
pansion will survive the functional integration, yields 
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The terms that survive the functional integral are of the 
form (ipoipo) n ipj 1 . . . tpj n tpj n - In the d — > oo limit, the 
large majority of these terms has j\ ^ j 2 ^ • • • ^ j n - 
There are (2n)\ of each terms of this type. Since each 
involves n pairs, and since there are d possible pairs to 
choose from, the total number of all such terms (note that 

j's are indistinguishable) is (2n)\(f) d ^°°> {2n)\d n /n\. 
Thus, in the large-c? limit ([5]) turns into 

'77S\2" d n (2n)\ 

(6) 
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where we have adopted the shorthand ip = tp(r) and ip 1 = 
ip(r'). This expression can now be inserted into Eq. ((4]), 
where the sum over n can be re-exponentiated, giving 



Z (1) (0) = / D(Vj>)e" So 



x exp 



d 2 i 



d 2 r'ipip' (ipip') 
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The superscript in Z^(0) signifies that this is the leading 
term in a large-c? expansion. Here we have defined fj = 
r\\fd as the new interlayer coupling, which remains finite 
as 7] — > and d — > oo. The j index has been dropped, 
since all layers are equivalent and are no longer coupled. 
The general correlation function is defined as 



(8) 



J D(ip,i/j)e- s o 
In the symmetric gauge, the correlation function in (J7J is 



(^(r)V(r')> 
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where z = (x + iy)/l is the complex coordinate within a 
single layer, I = ^J(Pq/2-kH is the magnetic length, and 
a is defined later. The integral in Eq. (J7|) is thus 
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The last equality follows from ip € LLL. 

Following Ref. we make change of variables 

tp(r) — ^JJ^z — Zi)e~' 2 ' / 4 , where {zi} are the positions 
of vortices. The interaction of {z^ is set by U^ 1 = 

y/{0A), where /3 A ({Zi}) = \ip\ 4 /\4>\ 2 is the Abrikosov 
ratio for arbitrary {zi} ( 7TT denotes a spatial average). 
The partition function for the zeroth layer becomes 



Z (1) (0) = / 

2tt1 2 sN 



S e fi = 



T 



dUe Ns We- s °* 



(11) 



a'W + 7^\®\ - Mn(27rPs|$| 2 ) 



Here N is the number of vortices {zi] and a' = a — fj 2 /a. 
The entropy function s(U) contains all the effects of lat- 
eral correlations among vortices {z^}, and knowledge of 
its exact form is equivalent to the exact solution for the 
thermodynamics of a single layer [l2j|. 

In the thermodynamic limit N — > oo, the saddle point 
method can be applied to integrals over $ and U in Eq. 
([TT]) . Minimizing with respect to |$| 2 gives 




(12) 



In order for this expression to be useful, we must deter- 
mine the form of a, as well as U. 

From Eq. ©, we have a" 1 = (2tt1 2 s/T) (-0(0)^(0)). 
Using this along with Eqs. (HJ and (TT2"1) . we obtain the 
following self-consistent expression for a: 



f3T 



a = a - — + (13) 

Solving this for a, and substituting the result into our 
expression for a' , we get 



1 - 



l + sgn(a),/l + 
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In solving for this expression, we must assume ft' = 
ft - 2fj 2 (Trl 2 sU 2 /T) > 0. ft' < leads to a < 0, which 
is clearly unphysical. The implications of ft' — > + at 
finite T are important and are discussed shortly. Eq. 
(TT41 constitutes our main theoretical result, allowing us 
to describe the system of coupled layers with a 2D GL- 
LLL action, albeit with a —¥ a' . Its innocent appearance 
notwithstanding, the change a — > a' actually entails an 
elaborate self-consistent calculation to determine the ul- 
timate dependence on T and H. Note that the next term 
in the large-ei expansion - arising from terms in ([5]) with 
one index repeated four times - modifies the quartic term 
ft in the 2D GL action. It is important to systematically 
incorporate such finite-ci corrections when addressing the 
details of interlayer correlations in real materials. 

Evaluating Eq. (fTTj) at its saddle point and using Eq. 
(TT21) . we obtain for the free energy density 
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where g = a' y / 2Trl 2 s/(2ftT) can be expressed using Eq. 
(IT!) as 
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In the above equation g a = \Js4>qH c2 (0) / 167tk 2 T c (0) , 

1 1/2 



S(h,t) 



1 



hi 



(t-t c2 {h)) 2 \2g 2 U 2 (g) 




and we used the GL result a 2 ,//? — H 2 2 (0)/8irn 2 [lj 

g(t, ft) (jTHJ) is the scaling variable of our theory. Since 
S(ft,i) depends on U(g), Eq. (TTBf has the form g = 
g(U(g)). U(g) is the same as in a purely 2D problem, 
but there g(t,h) — go(t — t c2 (ft.))/-\/fti, so the t and ft. 
dependencies in our case are very different. U(g) fol- 
lows from minimization of (jlll) and relies on knowledge 
of s(U). Here we can turn the problem around and ex- 
ploit the fact that /3 a (<?) interpolates between its high- 
and low-T limits of 2 and (3a = 1.159, respectively. In 
particular, 



U(g) =0.818-0.110 tanh 



+ C1 

C2 



(18) 



suggested in Ref. [12j, where c\ = 1.60 and c 2 = 2.66 from 
the fit to the Monte Carlo results of Ref. 

(H, yields a 

virtually exact solution for fluctuation thermodynamics 
20]. This expression for U(g) can then be used to solve 
self-consistently for g in Eq. (IPB"1) . 

It is now clear that the divergence in Eq. (|14l) . associ- 
ated with (3' 0+ and T -> T A = 2it1 2 t) 2 s/{P(3a), is en- 
dowed with special significance. As T is lowered toward 
Ta, g — > — 00 (since g oc a'), and thus J7(<?) — > \j\fB~A- 
Therefore, at /mi£e temperature Ta the system under- 
goes a Bose-Einstein condensation transition into the 
Abrikosov lattice state. In a purely 2D (fj — 0) the- 
ory, such a transition could occur only at T = 0. Once 
ij 7^ 0, this transition moves to finite Ta, which, over a 
large portion of an H — T phase diagram, is far below the 
vortex liquid-solid transition taking place at Tm, defined 

by g = g M 7 As H -s- 0, both T A and T M tend 

into T c2 (H). This echoes the phase diagram of layered 



21|. 



superconductors proposed in Ref. 

We now turn our attention to fluctuation thermody- 
namics [ijl, [23[. The magnetization follows from 
AttM = -(l/V)dF/dH, with |$ | 2 given in flT 



Att(/) sM 
T c (0) 



90 



(gU 2 - Uy/2 + g 2 uA . (19) 



Fig. [T] shows fluctuation magnetization data [9| for 
BaFe1.gCoo.2As2, and a fit of Eq. (IT9l) to the data. For 
this sample T c (0) = 23.6 K; and we obtain g = 5.8 us- 
ing the values H c2 (0) = 72 T for the upper critical field, 
k = 44 for the GL parameter 24 1, and s — 6.65 A for 
the interlayer spacing The demagnetization factor 

, which reduces the overall magnetization by a factor 
of 1 — Dm, is not known exactly for this sample, but can 
be estimated as Dm ~ 1 — nd/{2R), which is valid for a 
flat disk of radius R and thickness d R in a perpen- 
dicular magnetic field [25| . The sample used in Ref. 0] is 
rectangular in shape with length and width L ss 10d, so 
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FIG. 1. Scaled magnetization data from Ref. [9(, at fields 3, 
5, and 7 T, along with a theoretical fit from Eq. (|19|) . The 
theoretical scaling function ()19[) uses fitting parameters Dm = 
0.70 and fj/ao = 0.034 (solid line), with other parameters 
given in the text. The dashed line is the 2D case (77 = 0). 



we expect Dm ~ 1 — 7r/10. Fitting the data with respect 
to Dm and fj/a^, with other parameters fixed, yields the 
curve in Fig. [TJ and D M = 0.70, 77/^0 = 0.034. 

We now calculate the heat capacity C = -Td 2 F/dT 2 . 
From Eq. (fT5|) we obtain 



I dt 



hp. 



(2 3 [/ 2 -2t/v/2 + 5 2 C/ 2 ) 



1 

" 2t 




Here the heat capacity c = C/AC 2 d has been normal- 
ized to its 2D mean-field value, AC 2 d = Va^tKsPPA) = 



2VH c2 (0)g 2 t/((j>osf3 A ), and g is given by (TTB]). Fig. rj 
shows c for three different values of rj. As T — >• 0, there is 
a divergence in the specific heat, stemming from the fact 
that, for fj 7^ 0, g — 00 at /mite T — > Ta, as discussed 
before. This is suggestive of a nrst-ordcr Abrikosov tran- 
sition at Ta ; to describe its details our approach needs to 
be augmented either by the sixth order GL term (since 
j3' — > + at Ta) or finite d corrections, something left for 
future study. The specific heat, being a second deriva- 
tive, is rather sensitive to this divergence at low T, even 
for small fj, as we illustrate in the figure. 

Recent experiments on SmFeAsOi_ x F x (2(| suggest 
that the fluctuation conductivity follows an approximate 
2D scaling behavior of the form predicted by Ref. 14 1 
(see also Ref. [15J), where transport coefficients are de- 
rived from the time-dependent GL-LLL theory, within 
the Hartree-Fock approximation (/3a = 2). We follow 
Ref. 14j to obtain the fluctuation conductivity as 



Acr„ 



64A ob (0)£a6(0)T c (0)V sfo 



-£(<?), (21) 



where, in their case, the scaling variable g has its 



4 
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g (t + h-l)/Vht 

FIG. 2. Specific heat from Eq. J20J, with go = 3, h = 0.3, 
and t C 2(h) = 1 — h. The three curves have interlayer coupling 
values of fj/a = (solid), 0.002 (dashed), and 0.004 (dotted). 
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FIG. 3. Fluctuation conductivity data from Ref. [2y], along 
with a theoretical fit, with scaling variable g from (|16p . go = 
6.37, and fj/ao — 0.022 (solid line). The purely 2D curve (77 = 
0) is shown for comparison (dashed line). Other parameters 
are given in the text. 



2D form, i.e. 77 = in Eq. (p~5]) . In obtaining Eq. 
(f2"Tj) , we used the relations from the GL theory ao = 
e 2 H? 2 (0)\ 2 ab (0)/ m abC 2 K 2 and H c2 {0) = 0o/27r&(O), as 
well as the expression for the coefficient in the time- 
dependent GL equation Tq 1 ps 7rfia /8T c (0) [H[27j]. The 
scaling function in (|2ip has the form K.(g) = K,2D(g) = 
—g/2+^Jl + <? 2 /4, where now, of course, the scaling vari- 
able g must be changed to our Eq. (fL6|) . with 77 ^ 0. 

Comparison of the scaling function )C(g) to the experi- 
mental data in [26] is not straightforward since their sam- 
ple is a polycrystal. To compensate for this, we replace 
U(0) -> (U(0) 2 U0)) 1/3 , A ab (0) -> (A afe (0) 2 A c (0)) 1 /3 in 
the prefactor in (f2Tj) . Fig. [3] shows ZC(g) and the data for 
the optimally doped (x — 0.15, T c (0) = 51.5 K) sample at 
ff = 28 T. The coherence length is f o6 ( c )(0) = 24 (3) A 
2l; the penetration depth A ab(c) (0) = 2000 (16000) A 
28j; the upper critical field H c2 (0)/T c {0) = 7.8 T/K 
H fits snugly between \dH™ % /dT\ and \dH™ % /dT\ 
reported in Ref. 26]; and the interlayer separation 
s = 8.45 A 0. One' can see that the interlayer cou- 
pling leads to a strong enhancement of conductivity over 
its 2D form, even for modest values of fj/ao. 



In summary, we showed that a GL theory of coupled 
fluctuating superconducting layers in a magnetic field can 
be expressed as an effective, self- consistent single layer 
problem, in the limit of a large number of neighboring 
layers. Our approach can be generalized to other 2D, 
1+1D or 2+1D problems. Comparison of the theory with 
experimental results in the iron-pnictides is rather favor- 
able, and provides a means of making the quasi 3D nature 
of these materials more theoretically tractable. 
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